Internal Stress in a Model Elasto-Plastic Fluid 
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Plastic materials can carry memory of past mechanical treatment in the form of internal stress. 
We introduce a natural definition of the vorticity of internal stress in a simple two-dimensional 
model of elasto-plastic fluids, which generates the internal stress. We demonstrate how the internal 
stress is induced under external loading, and how the presence of the internal stress modifies the 
plastic behavior. 
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The internal stress, also known as the remanent stress 
in material science, is the stress which is maintained in 
a material without external mechanical supports. The 
internal stress is found in both discrete and continuous 
systems with widely different spatial scales. An examples 
of macroscale is the "tensegrity" of architecture pj > where 
the internal balance of the structural elements under ten- 
sion and those under compression maintains the solidity 
of the architecture. Its putative counterpart in biological 
cells has also been proposed 0, where the intracellular 
cytoskeletal networks are claimed to play the role of the 
structural elements. In lattice mechanics of solids, we 
could mention as examples of internal stress the lattice 
mismatch of crystal growth Q , the dislocations 0, |H El j 
the phase separation in alloys into the phases of differ- 
ent lattice constants 7\. In polymers, rubbers vulcanized 
under more than one deformed states (called the perma- 
nent set by Flory ||) would be an example of internal 
stress. The stresses in those systems, either in atomic or 
mesoscopic scales, are distributed and establish the me- 
chanical balance among them. The internal stress is thus 
a very common phenomenon. 

As the internal stress is not describable by the monova- 
lent and continuous elastic displacements Eshelby 
had introduced the "incompatibility tensor" [6j in the 
context of full three-dimensional lattice defects. Despite 
the several methods developed so far to describe the in- 
ternal stress, few studies has been done on the dynamics 
of the internal stress associated with plastic flows. Also 
little has been studied on the hysteresis related to the 
internal stress, especially on the plastic yielding under 
macroscopic deformations of the system including inter- 
nal stress 0, ^3 • I n this Letter we try to shed light on 
these problems through a simple model of elasto-plastic 
fluid in which we may use the analogue of the disloca- 
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FIG. 1: (a) System's geometry and the initial transient load- 
ing on the boundary (thick arrows) at x = 0. The fluid ve- 
locity du(x,y,t)/dt lies along the z-axis. (b) Time protocol 
of loading, s(t) (see the text). 



tion density 0,0 to characterize the internal stress. Our 
model is basically the Bingham model of plastic fluid 
supplemented with a finite shear compliance. This is a 
version of generalized Maxwell models of viscoelastic fluid 
with a longtime relaxation [1113. After describing our 
model, we will show how an internal stress is dynami- 
cally brought into a system under the external loading, 
and how the internal stress is maintained after the re- 
moval of the loading. Although the local threshold stress 
of plastic yielding (called the yield stress) is constant ev- 
erywhere, the internal stress maintained turns out to be 
generally subthreshold and inhomogeneous. This implies 
that the system with yield stress can bear the memory 
of past operations through the internal stress. Moreover, 
we will show that the presence of such internal stress 
modifies the plastic response of the system in a rather 
intriguing manner. 

Our model is constructed as follows. We constrain the 
fluid to flow in \\z direction while its flow profile may 
have heterogeneities in the lateral (x,y) directions. We 
denote the displacement of the fluid along the z-axis by 
u = u(x,y,t). The shear stress components are then 
essentially represented by the two-dimensional vector, 
°" = {^xz{x,y,t),a yz {x,y,i)). The evolution of u and a 
is given by the following isotropic continuum model: 
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where V = (d/dx,d/dy) and -F ext is the external force 
density, which we incorporate only as a boundary loading 
(see below), p and p are, respectively, the mass density 
and the elastic modulus per unit length along the z-axis. 
3>(<t) takes a two-dimensional vector value which is par- 
allel to er. This can be regarded as one of the generalized 
Maxwell models 0, > whose original form is recovered 
by replacing $(<r) by a. The special feature ascribed to 
the Bingham model [TlJ is that 3>(<r) includes a threshold 
(yield) stress, cry, of the plastic yielding: 



*(°0 



'0 (for < ov) 

<J — - — -<7y (for \<t\ > Oy) 



(3) 



The parameter r\ is the viscosity coefficient per unit 
length along the z-axis. Although not shown, this model 
can reproduce the one-dimensional propagation of the 
yielding (i.e. fluidization) front at the speed of the trans- 
verse elastic wave, v s — \J pj p. This process of fluid- 
ization, which is consistent with a numerical observa- 
tion [13, is not describable by the classical "Bingham 
model" where p = 00. 

Much more spectacular is the behavior in the higher 
dimensionalities, which can support the internal stress. 
In our model, we define the vorticity of the internal stress 
(VIS for short), which we denote by u>, as follows: 
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The VIS corresponds to the dislocation density in the 
continuum description of lattice defects |^, |5j , expressed 
in terms of stress. Our model @ gives the evolution of 
the VIS: 



ift = ~ divJ < 



(5) 



The first equation indicates that J is the conservative 
flux of u>. Eq. implies that uj is constant in time 
where there is no yielding, $ = 0. It also implies that 
the integral of VIS over a region is invariant if this is 
surrounded by a non-yielding region. 

To study the characteristics of the temporal evolution 
of the VIS as well as its influence upon the future plas- 
tic yielding, we have solved numerically Eqs. CQ~@- We 
will focus on the "hydrodynamic" or scaling regime whose 
spatiotemporal scales are much greater than the vis- 
coelastic timescale r = r\j ' p and the viscoelastic length- 
scale £ = v s t. To this end, we take the system of < x < 
L x = 120£ and \y\ < L y /2 = 50£, and the spatiotemporal 
mesh sizes, (Ax, Ay) = (0.390^,0.3910 and At = 5r/16. 
We apply a transient inhomogeneous stress o~ xz on the 
boundary at x — with a xz ^ for y ^ (Fig. ^a)), of 




FIG. 2: (a): The flux of the vorticity of the internal stress 
(VIS), J in units of cr Y /r, at t = lOOr. (b): The spatially 
distributed stress components (a xz ,a yz ) at t — 500r shown 
in the (a xz , <7 yz )-plane. The points on the curves give the 
stress components realized on the rectangular grids in the 
(x, j/)-plane, spaced by 15Ax and 16Ay, respectively. The 
dotted circle represents the yielding threshold, \cr\ = <jy- 
(c): The same data as (b) shown on the physical (x,y)- 
plane (in units of cty). (d): The VIS, u), for the same 
data as (b), shown by gray scales and contours. (Only 
the range of < x < 60£ is shown in (a,c,d), with the 
axes labeled in units of £.) The thick dots (•) in (c,d) 
at x = 17.0C, 26.8C (y > 0), 26. 8^ (y < 0), and 46.4£ indi- 
cate the locations where the subsequent yield starts under 
the specific uniform external shear stresses: (ct™* , °~^ff) — 
er ext (cos 6, sin 8) with 8 — n/2, 0, n, —n/2, respectively (see 
the text). 



the profile, a xz \ x=0 = fJ max s(t) tanh(y/L )- Here tr max = 
10 cry, Lo = 50 £, and s(t) is a slightly smoothened unit 
step-function, taking the value of unity only during the 
initial period of about 100t, with a smoothing period 
of approximately 1.5r (Fig. ^b)). The conditions on 
the other boundaries have been adjusted to mimic ap- 
proximately reflection-free walls 14]. For the temporal 
evolution we have used the two-step Lax-Wendroff (LW) 
scheme. The validity of this second-order scheme for the 
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present problem has been checked in the one-dimensional 
situation by comparing it with the Cubic-Interpolated 
Propagation (CIP) scheme of the third-order precision. 

Fig- a ) shows the snapshot of the flux of the VIS, J, 
at t = IOOt, the time at which the external loading is 
removed. The direction of J is governed by that of cr: 
from © and J5J), we have J _L er with \tr\ > ay. 

Figs.|3[b) and|2Ic) show the distribution of the inter- 
nal (or remanent) stress, cr = cr mt (x,y) = (a^(x,y), 
a^(x,y)), on the (a xz , cr yz )-plane and on the (x,y)- 
plane, respectively, at t — 500r, when all the plastic flow 
has died out and also the elastic waves have already left 
the system through the reflection-free boundaries (there- 
fore V-er lnt = 0). Evidently, the amplitude of the internal 
stress is by no means stuck to the threshold value, ay, 
but is distributed around zero. (The maximum value of 
er| in Fig.[2Ib,c) is about 0.515<ty-) This is because the 
stress field is elastically redistributed after the plastic de- 
formations so that it finally satisfies the static balance, 
V-er lnt = 0. In fact it is easy to show, from this condition, 
that the spatial average of er should be strictly equal to 
zero in the final state with no forces on the boundaries. 

Fig.EJd) shows the distribution of the VIS, uj (see (gj) 
calculated from the data of Fig.[2Ib,c). As compared with 
the spatial distribution of the internal stress, the VIS re- 
veals a peculiar localized peak, which is dislocated from 
the center of the noticeable large-scale circular distribu- 
tion of the stress itself |l5j . In the present setup of load- 
ing, the VIS (tu) is initially distributed rather broadly in 
the y-direction (not shown), reflecting the smooth spatial 
profile of the external loading at the boundary of x = 0. 
However, the flux J, which drives the VIS symmetrically 
with respect to the a;-axis (see Fig.^a)), acts to localize 
the VIS near this axis. The fan-shaped spatial charac- 
teristics of u> in Fig. E{d) suggests an approximate self- 
similarity in the distribution of the VIS. This self-similar 
nature can also be seen in the spatial Fourier spectrum 
of u) (Fig. I2J), which shows a power-law regime in the 
intermediate range of wave numbers. 

The state shown in Fig. E{b,c,d) retains the memory 
of the plastic flow in the past in the form of the VIS. 
This memory in turn influences the plastic response to 
the future external loading. As an example, if we apply 
quasi-statically a uniform external stress er — a CVLt t(9), 
where cr oxt > and t(Q) — (cos 6, sin 9), to the above 
mentioned state with an internal stress, the yield then 
occurs with a threshold amplitude of er ext which is gen- 
erally smaller than ay. Moreover, both the location (x, y) 
and the threshold amplitude of er oxt at which the yield 
starts depends on the orientation of the applied uni- 
form stress. As a demonstration, the thick dots in Fig. [21 
(c,d) indicate such locations for the external shear stress 
of several orientations. It is interesting to note that the 
locations where the yield starts is well distant from the 
region where VIS is concentrated. Moreover, the loca- 
tion of the first yield as a function of is not necessarily 
continuous nor single-valued, as seen by the two thick 
dots at x — 17.0 £ for the external stress: (0, er° xt ) with 
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FIG. 3: The log-log plot of the two-dimensional spatial 
Fourier spectral intensity of lo, integrated over the angle, vs 
the amplitude of the wave numbers (in units of go = (2£) _1 , 
the wave number beyond which the yielded fluid has no prop- 
agative modes). The exponent of the hydrodynamic regime 
(q < go) is read out to be about —1/3. 

ext > g 

We will describe how the system with the internal 
stress, cr lnt (x, y) (see Fig.|!Jb,c)), starts the plastic yield 
under further quasi-static application of a uniform exter- 
nal stress, er = a CKt t. First we note that, before the yield- 
ing, the uniform external stress is simply superposed onto 
the internal stress, <r mt (x,y), to give the actual stress as 
<r lnt (x,y) + er cxt i(6>). The external stress at the plastic 
yielding, which we denote by cry 1 , is, therefore, the solu- 
tion of the following equation: 

max \(T int (x, y) + afH(e)\ = ay, (6) 

where is fixed. Here those coordinates (x, y) that realize 
the maximal value on the left hand side tell the location 
(or locations) at which the plastic yielding starts. The 
threshold external stress thus obtained plays the role of 
apparent yield stress of the system with internal stress. 
Eq. Jfjjl is easily interpreted and solved graphically: As 
shown in Fig. E{b), the internal stress cr lnt (x,y) on the 
(a xz , (Tj /z )-plane forms a closed domain. The superposi- 
tion of an external uniform stress, <r cxt t((9), corresponds 
to a parallel displacement of this domain by a cxt t(9). 
If we increase tr ext , the displaced domain should eventu- 
ally touch from inside the circle of the yielding condition, 
er| = o~y (the dotted circle in Fig.|^b)). 

The above graphical solution immediately leads us to 
the following three qualitative conclusions: (i) The am- 
plitude of the apparent yield stress cr™' is generally less 
than ay, the value expected in the internal-stress- free 
system, (ii) If the apparent yield shear stress in one di- 
rection is less than ay, then it is true also in the opposite 
direction, to meet with the requirement of vanishing spa- 
tial average of the internal stress as mentioned above, 
(iii) The yielding point (x, y) which solves @l can be dis- 
continuous as the function of the direction of external 
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shear, 9. In fact, if we parallelly displace the domain 
of cr mt (x,y) by varying 8 while adjusting the value of 
tr cxt so that the domain remains tangent to the circle of 
|<x| = cry, the tangential point can undergo jumps. We 
expect a similar type of discontinuity to appear for the 
other criteria of plasticity as long as they share a similar 
mathematical structure to 10 ■ 

In this Letter we have defined and analyzed the inter- 
nal stress in a continuous model of elasto-plasticity, with 
the emphasis on the dynamics of the vorticity of the inter- 
nal stress ( VIS) as the source of the internal stress. Some 
of the insights obtained from our analysis could be gener- 
alized to other rheological models, and also to the models 
of higher dimensionality: (1) the internal stress after the 
yielding can be inhomogeneous, (2) the amplitude of the 
internal stress is generally less than the threshold value 
try, being redistributed so that the internal balance of 
elastic force is established, (3) the apparent threshold 
yield stress under subsequent global shear deformation 
is different from cry and anisotropic. In order to verify 
and develop these ideas, it is, firstly, essential to recog- 
nize the quantity like VIS which is responsible to the in- 
ternal stress and is invariant under elastic deformations. 
Secondly, it is important to analyze how such quantity 
is created and kept as the memory reflecting the past 
history of the system's plastic flow. Finally, we should 
analyze the consequences of the inhomogeneous internal 
stress such as the apparent threshold yield stress under 
the subsequent global loading. 

To be more realistic, the simple constitutive equation 



of the "Bingham type" should be modified so that it rep- 
resents the plastic yield as a slow dynamics involving the 
effect of temperature and the timescales of the observa- 
tion and the operation [tj 0, 0] . Since the definition 
of VIS is independent of the inertia effect, the concept 
of VIS should be applicable to other rheological mod- 
els in which the momentum can be ignored while some 
internal degrees of freedom of the material are relevant. 
Of particular interest would be the relation between in- 
ternal stress and slow relaxation: the (extended) distri- 
bution of relaxation time observed through the aging or 
creep experiments 01 could partly be attributed to the 
distribution of the internal stress. In systems consist- 
ing of discrete elements, like granular materials [T^j, a 
macroscopically homogeneous shear will be sufficient to 
cause internal stress among those elements, due to local 
anisotropy of elasto-plastic interactions. More quantita- 
tive analysis on this point could be made in future. As 
mentioned in the introductory part, the internal stress is 
found in many different systems not being limited to the 
rheological fluids. Systems composed of active elements 
(ex. biological cells), where the internal stress may be 
generated by the medium itself, would be highly of inter- 
est in this aspect. 
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